Deterministic ratchet from stationary light fields 
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Ratchets are dynamic systems where particle transport is induced by zero-average forces due to 
the interplay between nonlinearity and asymmetry. Generally, they rely on the effect of a strong 
external driving. We show that stationary optical lattices can be designed to generate particle flow 
in one direction while requiring neither noise nor driving. Such optical fields must be arranged to 
yield a combination of conservative (dipole) and nonconservative (radiation pressure) forces. Under 
strong friction all paths converge to a discrete set of limit periodic trajectories flowing in the same 
direction. 



The combination of broken space inversion symmetry 
and strong external driving is becoming an important 
principle for the controlled motion of a variety of ob- 
jects over a broad range of time and length scales. The 
ratchet effect [1, 2], whereby asymmetry and nonlinearity 
conspire to favor motion in one direction, has been real- 
ized in an increasing number of systems such as biological 
molecules [3], colloidal suspensions [4, 5], liquid droplets 
[6], superconducting vortices [7] and devices [8, 9], cold 
atoms [10], or magnetic domain walls [11]. A variety 
of questions has been explored theoretically, including 
the interplay between spatial asymmetries and thermal 
noise, time asymmetric forcing, topological driving, or 
collective phenomena [1, 2]. 

Despite the generality of these approaches, most pre- 
viously considered ratchets are based on some type of 
external driving, random or periodic. The underlying 
reason is that spatial asymmetry is not sufficient to in- 
duce directed motion: the necessary breaking of detailed 
balance requires also that the system be guided far from 
equilibrium. An external forcing is generally considered a 
necessary ingredient to produce the ratchet effect, which 
permits particle flow against friction in the absence of a 
net bias force. 

There is however a way of driving a system out of equi- 
librium which does not require external forcing and which 
has not been considered in the ratchet literature, namely, 
the use of nonconservative stationary force fields such as 
those resulting from radiation pressure vortices [12]. Op- 
tical fields are easily tunable in general and affect any 
polarizable object, from atoms to microscopic colloidal 
particles [13]. In this paper we show that a simple combi- 
nation of null-average conservative and nonconservative 
steady forces can rectify the flow of damped particles. 
A deterministic ratchet stemming from purely station- 
ary forces represents a novel concept in dynamics with 
considerable potential for fundamental and practical im- 
plications. For instance, it may permit the exploration 
of new forms of controlled atom motion in optical lat- 
tices [12, 14] and may be used to separate small particles 
with slightly different optical characteristics or friction 
properties [15, 16]. 



The use of radiation pressure forces is reminiscent of 
the drift ratchet, which relies on the effect of a surround- 
ing moving fluid [1]. However, drift ratchets have been 
based so far on the use of some external driving. Non- 
conservative steady forces have been contemplated in the 
context of chemical motors but only in cases where the 
molecule spatial coordinate moves in a conservative field 
[17, 18]. 

The two-dimensional field geometry arising in the in- 
tersection region of two standing waves oriented along 
the X and y axes illustrates a particular simple case of 
nonconservative optical forces. When the electric field is 
parallel to the z axis, the interference between the stand- 
ing waves pointing in perpendicular directions can lead to 
a lattice of radiation pressure vortices [12]. Specifically, 
for an electric field 



E(r) = Eq [sin(fca;) -|- i sin(fcy)] z , 



(1) 



the optical force on a polarizable particle becomes [12, 19] 



Fi = j£^o {V [sin^ikx) + sin2(fcy)] 
-I- 2/3V X [cos(/cx) cos(fc?/)z]} , 



(2) 



where a{Lo) — a' + ia" is the complex polarizability at 
frequency lo — ck (c is the speed of light in the surround- 
ing medium) and (3 = a" /a' . Notably, this force has a 
conservative (gradient) and a nonconservative (curl) con- 
tribution. The gradient term corresponds to the dipole 
force [20] , while the curl term accounts for radiation pres- 
sure [12, 21]. 

We consider the dynamics of a particle of mass m sur- 
rounded by a viscous medium creating a friction coeffi- 
cient 77. In the overdamped limit (V ma'EQk ^ 77) the 
equation of motion reads 



f = F(r)/77, 



(3) 



where F(r) is the stationary force field. We neglect here 
the fluctuating force due to thermal noise. This is jus- 
tified for low temperatures T such that ksT ^ I^I^'cj 
which is the scale of the energy barriers constraining the 
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particle motion. The general effect of a nonzero temper- 
ature will be that of softening the otherwise sharp sep- 
aration between the various flow regimes. Equation (3) 
allows us to treat the particle velocity at a given point 
as parallel to the local optical force. 

In Fig. la we show the contour lines for the potential 
generating the gradient force (first term) in Eq. (2) and 
the velocity flow map under the total force Fi given in 
(2). The square (-1, 3) x (-1, 3) is the unit cefl of the re- 
sulting periodic structure in the used reduced units. The 
vortex centers (where flux lines start) coincide with the 
potential maxima, i.e. the nodes of the electric field. For 
/3 < 1 pressure forces are too weak to free particles from 
their confinement within a quarter unit cell. In the case 
shown (/3 = 1.4), particles are able to reach the bound- 
aries of the quarter unit cells and move along the lines 
2kx/TT = 2n+ 1 and 2ky/TT = 2p+ 1 with n and p inte- 
gers. Due to the absence of activation barriers for /3 > 1, 
a minor degree of thermal noise suffices to permit parti- 
cles jumping between cells. As a result, particles diffuse 
throughout the lattice without a privileged direction. 

The long-time diffusion properties change qualitatively 
if a second laser is added within a suitable design. We 
consider the insertion of an independent laser which con- 
tributes a force 

F2 = d ^E^V [sm^k'x + <P)] . (4) 

The dimensionless parameter S gives the relative value 
between the intensity of the second laser and that of the 
conservative contributions of the first one [sec Eq. (2)]. 
Except for a final remark, we assume hereafter that the 
two lasers have the same wave number, k' = k. 

A most interesting case is that of (5 = 1 and phase shift 
(f) = 37r/8, which breaks the x —x inversion symmetry. 
Fig. lb reveals the dynamics under the total force field 
F = Fi + Fa = (Fr,Fy) for /3 = 1.4 [see Eqs. (2) and 
(4)]. The contour lines display the potential generating 
the gradient contribution to the total force F, while the 
flow map shows the total force (or velocity) field, includ- 
ing the nonconservative forces [curl term in Eq. (2)]. The 
vortex centers do not coincide with the potential maxima 
any more. Inspection of the fiow lines reveals the remark- 
able fact that all trajectories are open and point asymp- 
totically to the left: deterministic ratchet dynamics from 
stationary forces is possible. 

Figure 2 permits a more detailed analysis of the cir- 
culation patterns in the presence of the two lasers. The 
red/blue lines contain the points where the x/y compo- 
nents of the total force vanish. Their intersections define 
the fixed points of the flow map, where F = 0. Figures 
2a and 2d contain the same parameter values as Figs, 
la and lb, while Figs. 2b and 2c display intermediate 
regimes corresponding to (5 = 1 and /3 values 0.35 and 
0.45 respectively. 

A feature which is common to all maps in Figs. 1 and 2 



is that the vertical component of the total force vanishes 
along the horizonal lines 2ky/i: = 2p+ 1 with p integer, 
which prevents vertical diffusion. We are thus entitled 
to focus on the action in the interval 2ky/iT E (—1,3). 
On the contrary, the horizontal component vanishes at 
the vertical lines 2kx/'K = 2n -|- 1 (n integer) only in 
the absence of the second laser (S = 0, sec Figs, la and 
2a). The most important effect of the second laser is that 
of rendering nonzero at the vertical lines separating 
the various quarter unit cells. This is clear in Figs, lb 
and 2b-d. Although the particle polarizability is not the 
easiest parameter to tune in a typical experiment, in the 
following we study the circulation trends as a function of 
P because, against the background of a flxed conservative 
force, a clearer analysis is possible. 

Figure 2b shows a case where nonconservative forces 
are too weak (/3 = 0.35) to permit diffusion. Parti- 
cles can move across squares but they always end up 
in an attractive equilibrium point. There are two types 
of them, located slightly to the right/left of represen- 
tative points P/Q in the figure. An important change 
occurs for j3 = 0.393, which is the root of the equation 
F^{x,tt/2) = 0, d^F^{x,Tr/2) = 0. At that value of /3 
the system undergoes a saddle-node bifurcation, as the 
attractive fixed points near P merge into point P and 
disappear; the short horizontal stretch comprising point 
P, with velocity pointing to the right, shrinks to zero. 
However, for values of /? only moderately above this first 
critical value (0.393), diffusion is still not possible. This 
is clearly shown in Fig. 2c (/3 = 0.45), where all trajecto- 
ries die in the attractive fixed points of the near-Q type, 
which still remain. 

A major, qualitative change occurs at the critical value 
/?c = 1.03, which is the root of the equation Fx{x, 7r/2) = 
0, dyFy{x, Tr/2) = 0. At /? = /3c there is a subcritical 
pitchfork bifurcation changing the stability of the repre- 
sentative fixed point near Q, which is not attractive any 
more but rather becomes a saddle point (attractive in the 
X direction and repulsive in the y direction). A typical 
case is shown in Fig. 2d, where (3 = lA. The result is 
that no attractive fixed points exist any longer. Particles 
can move throughout the lattice but, very importantly, 
only towards the left and in a deterministic manner (i.e. 
without requiring thermal activation). From stationary 
forces of null average we have designed a ratchet system 
which requires neither noise nor driving. 

The deterministic character of the resulting ratchet is 
best appreciated in Fig. 3, where 100 trajectories are 
shown with initial positions chosen at random within the 
unit cell. From a generalization to two-dimensional man- 
ifolds of the Poincare-Bendixson theorem [22] , and as can 
be clearly seen in Figs. 2d and 3, limit attracting periodic 
trajectories exist. These trajectories satisfy the differen- 
tial equation dy/dx = Fy/F^ with periodic boundary 
conditions. 

Figure 4a shows the average velocity of the limit trajec- 
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tories as a function of the parameter /?, normalized to the 
amphtude of the total force (proportional to 2/3 + 1). The 
minimum value /3c = 1-03 is required to have a nonzero 
mean velocity (ratchet effect). A detailed analysis shows 
that, for /3 slightly above /3c, the average velocity grows 
linearly, v{(3) ^ P — (3c. For large values of /3, the velocity 
grows slower than hnear, v ~ 7r/3/21n/3. 

Figure 4b shows the average velocity for fixed /3 = 1.4 
as a function of the parameter S, which gives the relative 
strength of the second laser and is thus easier to tune. 
The limiting behaviors are v ~ {\nS)~^ for S 0+, and 
V ~ for S ^ 5~, where ^c = 3.29 gives the maxi- 

mum intensity of the second laser permitting rectification 
of motion [see Eq. (4)]. For S > Sc conservative forces 
dominate and the ratchet effect disappears. 

Wc have focused our study on the case k' = k and 
(5 = 1 in Eq. (4), which has the advantage of involving a 
single optical wavelength. Other conservative forces are 
of course possible. We note here that the case k' = 2k, 
5 = 1/4, and = 0, which often stands as a paradigmatic 
ratchet potential [1] , also yields a deterministic rectifica- 
tion of motion under the action of radiation pressure from 
the optical vortex lattice. 

We conclude that suitably designed optical lattices can 
provide a playground for novel forms of nonconserva- 
tive stationary dynamics. As a remarkable example, wc 
have demonstrated theoretically that deterministic ratch- 
ets from stationary light forces are possible. Figure 3 of- 
fers a vivid picture of how properly tailored steady fields 
can drag particles without a net bias force. The fine 
sensitivity of the flow regime to polarization and fric- 
tion properties can pave the way to novel protocols of 
separation of different particle species [15, 16]. We may 
generally expect the new dynamic concept here unveiled 
to open new vistas in the controlled transport of polariz- 
able neutral bodies such as cold atoms or small dielectric 
particles. 
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FIG. 1: Equipotential lines are plotted with interstitial col- 
ors growing from blue (potential minima) to white (maxima) . 
Red arrows show the individual trajectories of particles sub- 
ject to the total force field, which includes the conservative 
force deriving from the shown potential and a nonconservative 
force obtained from the curl of a not shown vector potential. 
The square (—1, 3) x (—1, 3) is a unit cell of the resulting pe- 
riodic structure. Points P and Q are described in the caption 
of Fig. 2. (a) Single-laser case [5 = 0,/3 = 1.4 in Eqs. (2) 
and (4)]; particles are confined within a quarter unit cell but 
even a minor amount of thermal noise allows them to diffuse 
isotropically throughout the periodic structure, (b) Double- 
laser (5 = l,/3 = 1.4); all trajectories evolve to the left, show- 
ing a ratchet effect which stems from purely deterministic and 
stationary forces. 
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FIG. 2: Grey arrows show particle trajectories under the effect 
of the total (nonconservative plus conservative) force. The 
force is horizontal/ vertical at the blue/red lines, (a) 5 = 0, 
/3 = 1.4; (b) (5 = 1, /? = 0.35; (c) 5 = 1, /? = 0.45; (d) 5 = 1, 
P — 1.4. Points P and Q become critical (zero force) at the 
transition between flow regimes, defined by the presence or 
absence of the various classes of attractive fixed points. 
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FIG. 3: A set of 100 trajectories is shown for the setup of 
Figs, lb and 2d. Initial positions are random within the 
chosen unit cell. AU trajectories converge to a limit periodic 
trajectory which flows to the left. 




FIG. 4: Average velocity v of the limit periodic trajectories, 
in units of o'Eik/Arj [see Eqs. (2)-(4)]. (a) v/{l + 2/3) as 
a function of /3 for second- laser intensity 5 — 1. (b) w as a 
function of 5 for /3 = 1.4. 



